Cosmological anisotropy from non-comoving dark matter and dark energy 
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We consider a cosmological model in which the two major fluid components of the Universe, 
dark energy and dark matter, flow with distinct four-velocities. This cosmological configuration is 
equivalent to a single anisotropic fluid, expanding with a four-velocity that is an appropriate com- 
bination of the two fluid four- velocities. The energy density of the single cosmological fluid is larger 
than the sum of the energy densities of the two perfect fluids, i.e., dark energy and dark matter, 
respectively, and contains a correction term due to the anisotropy generated by the differences in the 
four-velocities. Furthermore, the gravitational field equations of the two-fluid anisotropic cosmolog- 
ical model are obtained for a Bianchi type I geometry. By assuming that the non-comoving motion 
of the dark energy and dark matter induces small perturbations in the homogeneous and isotropic 
Friedmann-Lemaitre- Robertson-Walker type cosmological background, and that the anisotropy pa- 
rameter is small, the equations of the cosmological perturbations due to the non-comoving nature of 
the two major components are obtained. The time evolution of the metric perturbations is explicitly 
obtained for the cases of the exponential and power law background cosmological expansion. The 
imprints of a non-comoving dark energy - dark matter on the Cosmic Microwave Background and on 
the luminosity distance are briefly discussed, and the temperature anisotropies and the quadrupole 
are explicitly obtained in terms of the metric perturbations of the flat background metric. There- 
fore, if there is a slight difference between the four-velocities of the dark energy and dark matter, 
the Universe would acquire some anisotropic characteristics, and its geometry will deviate from the 
standard FLRW one. In fact, the recent Planck results show that the presence of an intrinsic large 
scale anisotropy in the Universe cannot be excluded a priori, so that the model presented in this 
work can be considered as a plausible and viable working hypothesis. 

PACS numbers: 04.50.Kd, 04.20.Cv, 04.20.Fy 



I. INTRODUCTION 



The recently released Planck data of the 2.7 full sky 
survey [l|-|5[ have shown a number of intriguing features, 
whose explanation will certainly require a deep change 
in our standard view of the Universe. These recent ob- 
servations have measured the cosmic microwave back- 
ground (CMB) to an unprecedented precision. One of 
the most striking results of the nominal mission is that 
the best-fit Hubble constant has the value Hq = 67.4±1.2 
km s^Mpc -1 , with a dark energy density parameter 
SIa = 0.686 ± 0.020, and a matter density parameter 
Q, M = 0.307 ± 0.019. Even though generally the Planck 
data confirm the foundations of the ACDM (ACold Dark 
Matter) model, the observational data show some ten- 
sion between the fundamental principle of the model and 
observations. For example, the Planck data combined 
with the WMAP polarization data show that the index 
of the power spectrum is given by n s = 0.9603 ± 0.0073 
@, Q: at the pivot scale fco = 0.05 Mpc -1 , which rules 
out the exact scale-invariance (n s = 1) at more than ha. 
In addition to this, the joint constraints on r (tensor-to- 
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scalar ratio) and n s can constrain inflation models. For 
example, inflationary models with a power-law potential 
</> 4 cannot provide reasonable number of e-folds (between 
50-60) in the restricted space of r — n s at around a 2c 
level. 

More interestingly, the search for the background ge- 
ometry and topology of the Universe reveals that a 
Bianchi pattern, corresponding to an isotropic geometry 
of the Universe may account, in a quite efficient way, for 
some large-scale anomalies seen in the Planck data Q . A 
Bayesian search for an anisotropic Bianchi Vllh geometry 
was performed, by using the Planck data in [5j . In a non- 
physical setting, where the Bianchi parameters are de- 
coupled from the standard cosmology, the observational 
data favour a Bianchi component with a Bayes factor of 
at least 1.5 units of log-evidence. However, in the physi- 
cally motivated setting where the Bianchi parameters are 
fitted simultaneously with the standard cosmological pa- 
rameters, no evidence for a Bianchi Vllh cosmology was 
found @. 

According to the ACDM model, the material content 
of the Universe consists of two components: dark energy, 
and dark matter, respectively jg]. The observed late-time 
cosmic acceleration of the Universe [7j can be successfully 
explained by introducing either a fundamental cosmolog- 
ical constant A, which would represent an intrinsic cur- 
vature of space-time, or a dark energy, a hypothetical 



fluid component in the form of a zero-point-energy that 
pervades the whole Universe, which would mimic a cos- 
mological constant (at least during the late stage of the 
cosmological evolution) . Currently one of the main dark 
energy scenarios is based on the so-called quintessence 
|8|, where dark energy corresponds to a scalar particle 
<t>. There are also other dark energy candidates, such as 
the Chaplygin gas model [9(, in which dark energy is a 
hypothetical substance that satisfies an exotic equation 
of state in the form p = A/ p~ a , where p is the pressure, 
p is the energy density, and A and a are positive con- 
stants. Modifying gravity at large scales, or considering 
extra-dimensional cosmological models can also provide 
an explanation of the late acceleration of the Universe 

Mm 

Dark matter is thought to be composed of cold neu- 
tral weakly interacting massive particles, beyond those 
existing in the Standard Model of Particle Physics, and 
not yet detected in accelerators or in direct and indirect 
searches. There are many possible candidates for dark 
matter, the most popular ones being the axions and the 
weakly interacting massive particles (WIMP) (for a re- 
view of the particle physics aspects of dark matter see 
[12j]). Their interaction cross section with normal bary- 
onic matter, while extremely small, are expected to be 
non-zero and we may expect to detect them directly. 
Scalar fields or other long range coherent fields coupled 
to gravity, have been considered as potential dark matter 
candidates [lSq . Recently, the possibility that the galac- 
tic dynamics of massive test particles may be understood 
in the context of modified theories of gravity, without the 
need for dark matter, was also explored [13 |. 

Now, precise observations of the cosmic microwave 
background radiation enable us to test the fundamental 
predictions of inflation on primordial fluctuations such 
as scale independence and Gaussianity QHa]. The sta- 
tistical isotropy has been a robust prediction protected 
by the cosmic no-hair conjecture, which claims that in- 
flation washes out classical anisotropy. However, several 
recent observations of the large scale structure of the Uni- 
verse have questioned the principles of homogeneity and 
isotropy [ljj . As combined with the recent Planck re- 
sults, these robust observational results show that the 
presence of an intrinsic large scale anisotropy in the Uni- 
verse cannot be excluded a priori. Recently, the pos- 
sibility that the geometry of the Universe may not be 
of the Friedmann-Lemaitre-Robertson- Walker (FLRW) 
typ e was investigated from different points of view 
[la . |l6j. In particular, in [15| the homogeneous locally 
rotationally symmetric (LRS) class of metrics were con- 
sidered as describing the global geometry of the Universe. 
Such spacetimes possess a preferred direction in the sky 
and at the same time a CMB which is isotropic at the 
background level. A distortion of the luminosity dis- 
tances has been obtained, and used to test the model 
against the CMB and supernovae data. It has also been 
shown that the latter exhibit a higher-than-expected de- 
pendence on angular position. Local constraints from 



limits from simulated distance measurements for various 
distributions of survey fields in a Bianchi I anisotropic 
universe were obtained [lTj • The effects of fitting for line 
of sight properties where isotropic dynamics is assumed 
were also considered, as well as the comparison of sen- 
sitivities of observational probes for anisotropies, from 
astrophysical systematics and from dark energy. 

In addition to these, especially in the view of the latest 
Planck data, the empirical evidence for the existence of 
cosmological anisotropies is strong. But, its physical ori- 
gin is still unknown, although the most popular explana- 
tion being that the primordial fluctuations are deviating 
from isotropy [ij. However, there is no clearly established 
physical mechanism that could lead to such deviations. 

It is the purpose of the present paper to consider 
another possibility of explaining the presence of the 
anisotropies on a large cosmological scale, namely, the 
possibility that the two major components of the Uni- 
verse, dark energy and dark matter, are not comoving. In 
standard cosmology the assumption that all constituents 
of the Universe move in the same frame, with the same 
four- velocity, is a fundamental (but not rigorously justifi- 
able) assumption. This allows adopting a frame, which is 
comoving with matter in all its forms, in which the com- 
ponents of the four velocities u M of all components can 
be chosen as u M = (1,0,0,0). Consequently, the global 
thermodynamic parameters of the Universe are just the 
sum of the individual thermodynamic parameters. Thus, 
for example, in the standard model the energy density 
of the Universe p is given by p = X)"=iPjj where pi, 
i = 1,2, ..,n are the energy densities of the individual 
components. 

However, there is no basic physical or observational 
principle that would require that all components in the 
Universe should have the same four-velocity, that is, 
one could describe the dynamics of the Universe in a 
single frame, comoving with all the constituent compo- 
nents. In this paper we will investigate the situation 
in which the Universe may be modeled as a mixture of 
two non-interacting perfect fluids, namely dark energy 
and dark matter, respectively, possessing different four- 
velocities. This configuration is formally equivalent to 
a single anisotropic fluid |18l - l21| . Therefore, if there is 
a slight difference between the four- velocities of the dark 
energy and dark matter, the Universe would acquire some 
anisotropic characteristics, and its geometry will deviate 
from the standard FLRW one. 

The thermodynamic parameters (energy density and 
pressure) of the non-comoving two-fluid Universe are ob- 
tained by reducing the system to a single fluid, described 
by an anisotropic energy-momentum tensor. The energy 
density of the single fluid is greater than the sum of the 
energy densities of the dark energy and dark matter, re- 
spectively, and it contains a correction due to anisotropy. 
As a cosmological application of the model we consider 
an effective induced Bianchi type I geometry, generated 
by the non-comoving expansion of the dark energy and 
dark matter, respectively. The general properties of the 



system are discussed, and the expression of some obser- 
vationally important parameters (expansion and shear) 
are obtained. By using a simple perturbative approach, 
some simple cases of non-comoving dark anergy and dark 
matter cosmological models are constructed. 

The possibility that there may exist different rest 
frames for dark matter and dark energy has been also 
investigated in [22f , where it is shown that the existence 
of large scale bulk flows could be signaling the presence of 
moving dark energy at the time when photons decoupled 
from matter. Bianchi I type cosmological models, the 
evolution of the homogeneous perturbations up to sec- 
ond order, and the effects on the CMB quadrupole were 
also considered. 

The possibilities of observationally testing the non- 
comoving dark energy-dark matter are briefly discussed. 
In particular we analyze the imprints of the anisotropy 
on the CMB spectrum, and on the luminosity distance. 
The expressions of the distribution of the temperature 
anisotropies and of the quadrupole C2 are explicitly ob- 
tained in terms of the metric perturbations of the back- 
ground flat geometry. This opens the possibility of ob- 
servationally constraining the parameters of the model 
by using the recent Planck results [lj , as well as the type 
la supernovae data Q . 

The present paper is organized as follows. The reduc- 
tion of a two-fluid system to a single effective anisotropic 
fluid description is discussed in Section [TTJ The cosmo- 
logical gravitational field equations of the model are ob- 
tained in Section irm where some general properties of the 
non-comoving two fluid Universe are also discussed. The 
evolution of the cosmological perturbations due to the 
non-comoving motion of the dark energy and dark mat- 
ter are considered in Section IIVI The observational im- 
plications of our results are investigated in Section [V] We 
discuss and conclude our results in Section [VlJ Through- 
out this paper, we use the natural system of units with 
c = 8ttG = 1. 



II. ANISOTROPIC FLUID DESCRIPTION OF 

THE NON-COMOVING DARK ENERGY-DARK 

MATTER COSMOLOGICAL MODELS 

We assume that the Universe consists of a mixture of 
two basic fluid components: dark energy, with energy 



density and pressure pde and pde, and four-velocity 
respectively; and dark matter, with thermody- 



L DE> 



namic parameters pdm and pdm, respectively, and four- 
velocity u DM . In standard cosmology it is assumed that 
the two cosmological fluids are comoving, which implies 



L DE 



L DM 



(1) 



This condition allows to choose a comoving frame in the 
study of the cosmological dynamics, where the compo- 
nents of the four velocity of all components can be re- 
duced to the form m m = (1,0,0,0). However, in the 
present study we relax this condition, by assuming that, 
at least in some periods of the cosmological evolution, the 
two fundamental components of the Universe may have 
had different four-velocities, so that 



"DE 



r u dm- 



(2) 



The physical and thermodynamical properties of the 
two-component cosmological fluid arc described by the 
energy-momentum tensor T^ u of the Universe, given by 
the sum of each individual components as 



T^ v = (pde + Pde) u de u de - p D E9^ v + 

(Pdm + Pdm) u dm u dm - pdm9^ ' ■ (3) 



The four-velocities of the dark energy and dark mat- 
ter arc normalized according to u de ude ^ = 1 and 
u dm Udm m = 1' respectively. In the following, we as- 
sume that the energy density and the pressure of the 
dark energy satisfy the conditions pde + Pde 7^ 0, that 
is, dark energy cannot be described by a "pure" cosmo- 
logical constant. 

The study of the cosmological models described by 
an energy-momentum tensor having the form given by 
Eq. ([3]) can be significantly simplified if we cast it into 
the standard form of perfect anisotropic fluids. This can 
be done by means of the transformations y/^ 
and u DM 



T>Af> 



respectively, so that [18|-|20( 
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'■DE _ 
L DM 



cos a 



Pde+Pde 
Pdm+Pdm 



sin a 



Pdm+Pdm 
Pde+Pde 



cos a 



sma 



'DE 
l DM 



(4) 



represents a "rotation" of the four- vector velocities in mations (j4j take the form 



the {u DE ,' 



■DM 



) velocity space. Explicitly, the transfor- 



ms 



h DE 



L DE' 



Pdm + Pdm . 
Pde + Pde 



(5) 



L DM 



h DM 



L DM ' 



PDE 



Pdm + Pdm 



PDE a 

u EE sm a . 



(6) 
respectively. 

The velocity transformations given by Eq. (j4]) 
leave the quadratic form (pde + Pde)u ee u ee + 
(Pdm +Pdm)u d . m u v dm invariant. Thus, 



T ia> ( U P 



DEt u DM, 



= T^ (i 



■DEi U DMJ 



(7) 



Next, we choose u*p E and u* EM such that one becomes 
timelike, while the other is spacelike. Therefore, the 
two transformed four-vector velocities satisfy the condi- 
0. With the use of Eqs. ©-© and 
0, we obtain first the rotation angle as 



tion u*g E vr DM 



tan 2a 



\J(pde +Pde) (pdm +Pdm) 
Pde + Pde ~ (pdm + Pdm) 



'de u dm , 



(8) 

Note that if pde + Pde = 0, which corresponds to the 
case of the cosmological constant, the rotation angle is 
a = 0, and therefore the rotation in the velocity space 
reduces to the identical transformation. 
Now, by defining the quantities 



V* = 



X" = 



"DE 



V 



U DE U a DE 



"DM 



v- 



(9) 

(10) 



"DM "-DM a 



2 (pde + Pdm - Pde - Pdm) + 



and 



T^V^V V 



(Pde + Pde) u*£ e u* de a - (p DE +Pdm) , (11) 



= (pde + Pdm) - (pdm + Pdm) u*d M W^m J, 12 ) 



n = pde +pdm, 



(13) 



respectively, the energy-momentum tensor of the two 
non-interacting perfect fluids can be written as 

t^ = (e + n) y"V" - n 9 ' w + (* - n) x ^x v , (14) 

where V% = 1 = -X^X^ and x^% = [HI3. 
The energy-momentum tensor of the two-fluid cosmo- 
logical model, in which the components have different 
four- velocities, given by Eq. (ThO) . is the standard form 
for anisotropic fluids 20]. 

The energy density e and the radial pressure ^> of the 
anisotropic Universe are given by 



~\l (pDE +PDE + PDM +PDAl) +i(pDE + Pde) (PDM + Pdm) 



(u ee Udm h) — 1 



(15) 



and 



* 



(pde + Pdm - Pde - Pdm) + 



~ V (pde +Pde - Pdm -Pdm) 2 +^(pde +Pde) (pdm +Pdm) (u^udmji) , 



(16) 



respectively [18H20J. 



gives first 



Since the two cosmological fluids, the dark energy, and 
the dark matter, have different four-velocities, we may 
write u^ e udm « = 1 + b/2, where generally b is an ar- 
bitrary function of the radial and time coordinates. The 
functional form of b can be obtained from Eq. ([5J , which 



u EE u DM ^ = -tan 2a 



Pde +Pde 
Pdm +Pdm 



PDM + PDM 

Pde + Pde 



while b can be obtained as 
b = 



Pde + Pde 



PDM + PDM 



Pdm + Pdm 



Pde + pde 



(17) 



tan2a-2. (18) 



By assuming that the physical parameters of the dark 
energy and dark matter satisfy the condition 



26 1 



(Pde + Pde) (pdm + Pdm) 
8 / (pde + Pdm + Pde + Pdm) 



« 1, (19) 



the energy density, the radial and tangential pressures of 
the Universe can be obtained as 



Pde+Pdm- 



* 



Pde+Pdm- 



b 

b 

F 2 



1 



(P 



DE 



Pde) (pdm +Pdm) 



(pde + Pdm + Pde 
(Pde + Pde) (pdm + Pdm) 



Pdm) 
(20) 



(P 



de + Pdm + Pde 



and 



II = PDE +PDM, 



Pdm) 

(21) 

(22) 



respectively. As one can readily verify from Eq. (I20|) , the 
energy density of the non-comoving two fluid Universe 
is different from the sum of the energy densities of the 
component fluids. 



III. ANISOTROPIC TWO FLUID 
COSMOLOGICAL DYNAMICS 

In the present Section, we consider some cosmological 
applications of the non-comoving dark energy-dark mat- 
ter cosmological model. As a first step in our analysis, we 
write down the gravitational field equations correspond- 
ing to the anisotropic evolution of the Universe. The 
general properties of the model are analyzed, and it is 
shown that the non-comoving nature of the cosmologi- 
cal dynamics affects the observational parameters of the 
Universe, by inducing some specific anisotropic effects in 
the evolution of both expansion and shear parameters. 
Finally, an explicit dark energy model is considered, and 
the general solution of the field equations, describing the 
dark energy-dark matter mixture is obtained. 



where W v = g^ — V^V U . The conservation of the 
energy-momentum tensor V M T^ = yields 

e + (e + II) Vj.V" - (tf - n) x ,fi V„ = 0, (26) 



(e + n) v" + /i^v^n 
+ (* - n) (v„ x "x M - 



(tf-n)V 
•x>^) = o, 



(27) 



where the overdot and the prime are defined as () = 
V M () V", and ()' = V M () x M , respectively. 

Since the two-fluid dark energy-dark matter cosmo- 
logical mixture represents an anisotropic fluid, on a 
cosmological scale the corresponding geometry is also 
anisotropic. In comoving Cartesian coordinates x° = i, 
x 1 = x, x 2 = y, and x — z we may choose V 1 = 
V 2 = V 3 = 0, V°V = 1 and X ° = X 1 = X 2 = 0, 
X 3 X3 = — 1) respectively }18l420| . Therefore the com- 
ponents of the energy-momentum of two non-comoving 
cosmological perfect fluids take the form 



J- n 



: £, 



Tl = T 2 



-n. 



T$ 



-*. 



(28) 



where e is the total energy-density of the mixture of 
fluids, $ = P z is the pressure along the z direction, 
while II = P x — P y is the pressure along the x and y 
directions. The energy-momentum tensor components 
given by Eqs. (|2"5)l shows that the pressure ^ along the 
z— direction may be different with respect to the pressure 
n on the x — y plane. 

In a homogeneous Universe with a two-fluid mixture 
the simplest geometry exhibiting this symmetry of the 
energy-momentum tensor is the flat Bianchi type I ge- 
ometry, with the line element given by 



ds 2 = dt 2 - a\(t)dx 2 - al(t)dy 2 - al(t)dz 2 



(29) 



where ai(t) (i = 1,2,3) are the directional scale factors. 
In this geometry the gravitational field equations take 
the form 



'AH + Hi + Hi 



Hi 



1 



(e + * + 211) ; 



(30) 



A. Gravitational field equations of the two 
non-comoving fluid cosmological model 

For the two fluid energy-momentum tensor, with com- 
ponents given by Eqs. (fH)l . the Einstein gravitational 
field equations can be written as }19l | 



RauV^V V 



(e + 2IT + ip) . 



(23) 



and 



^(VF i ) = i(e-n), a = l,2, (31) 



viw-ki'-v* 



(32) 



respectively. We have denoted 



^ifj.i/hg.h^ 



RarVhZ = 0, 



e - - (211 + tt) 



I a A 



(* - n) ( XaXX - g^A ) , (25) 



(24) 



V = a 1 a 2 a 3 , H t = — , 



and 



z = 1,2,3 



3 \ 2-^, % \ 31/ 



(33) 



(34) 



respectively. The expansion parameter 9 is given by 
9 = 3H. From the conservation of the total energy- 
momentum tensor of the two-fluid cosmological model 
we obtain the evolution equation of the energy density of 
the system, given by 



£ + 3 (e + n) H - (* - n) H 3 = 0. 



(35) 



General properties of the anisotropic Bianchi I 
two fluid cosmological model 



B. 



From Eqs. (|31[) it follows that for the considered form 
of the anisotropic energy-momentum tensor one can take 
ai(i) = 0,2 (£) without any loss of generality. Therefore 
the expansion parameter of the Universe can be obtained 

as 

9 = 2^ + ^, (36) 

a\ a 3 

while the shear scalar a of the Universe is given by 



1 



V3 \ a 3 



«3 



Ol 

a i 



Equations (1361) and (|37j) allow the determination of the 
directional Hubble parameters -Hi and H3 in terms of 
observable cosmological parameters as 



Hi= 3-7r 



^=3 + 7r' (38) 



respectively 

The large-time behavior of an anisotropic cosmological 
model can be obtained from the study of the anisotropy 
parameter A, defined by 



\Y. 



; 'Hi-H^ 2 



3^ V H 

l—l 



(39) 



If A = the cosmological model is isotropic. The 
anisotropic cosmological pressure differences can be ob- 
tained, in terms of the observable quantities, as 



* - II = VSa + 9a, 



(40) 



(37) or, equivalently, 



Via + 9a 



(pDE + RDM) + A/ {PDE + PDE - RDM - PDAl) + 4 (pDE + PDe) (PDM + PDAl) (Up E U D M fi) 



r 



ill) 



By assuming that the parameter /3(i), describing the 
deviations from a comoving motion, is small, then we 
obtain 



(pde +Pde) (pdm +Pdm) (up E u DM M ) 
(pde + Pde + Pdm + Pdm) 



(PnE + Pde + Pdm +Pdm) 

(42) 

To first order in 0, we obtain the following expressions 
of the energy density and z-direction pressure of the two 
non-comoving fluid cosmological model, 

£ « Pde + Pdm + P (pde + Pde + Pdm + Pdm) , (43) 

and 

* = Pde + Pdm + P {Pde + Pde + Pdm + Pdm) , (44) 

respectively. Therefore the total energy density and pres- 
sure of the two fluid cosmological model is not equal to 
the sum of the energy densities and pressures of the in- 
dividual matter components of the system. 
By adding Eqs. (|31[) and (|32|) we obtain 



H+3ff2 4 £ 4( n+ ! 



(45) 



or, equivalently, 
V- 



n+i*-| e |r = o. 



(46) 



IV. PERTURBATIVE APPROACH TO 

NON-COMOVING TWO FLUID 

COSMOLOGICAL MODELS 

In the following, we consider a simple perturbative 
solution of the field equations (|30p - (f3"2")l by assuming 
that the anisotropic expansion due to the non-comoving 
motion of the dark energy and dark matter represents 
a small perturbation of the background isotropic, flat 
FLRW geometry, described by the isotropic scale factor 
a(t). We also assume that the expansion rates along the 
x and y axes are the same, so that a\(t) — ai{t), and the 
perturbations along these axes are identical. Therefore 
the scale factor of the perturbed Bianchi type I geometry 
can be given as 



(kit) = a(t) + 6i(t), 



1,2,3, 



(47) 



where 5i(t), i — 1,2, 3, are small corrections in the scale 
factor a(t) of the isotropic Universe, induced by the small 



differences between the four- velocities of the dark matter 
and dark energy. Moreover, we assume that 



5i(t) = 6 2 (t). 



(48) 



The Hubble parameter of the background isotropic space 
is H (t) — a/ a. For the directional Hubble parameters 
of the anisotropic expansion we obtain 



_ a(t)+6i(t) ^ 
a(t) + di(t) 



1 



Si(t) 



a(t) 



Sj(t) 

ait)' 



(49) 



where i = 1,2,3. Thus, the mean Hubble factor can be 
obtained as 



H(t) = l 



3tf (t)-F (t)M + M 

a(t) a(t) 



where we have denoted 

S(t) = 26 x (t) + 6 3 (t), 



(50) 



(51) 



while the square of the directional Hubble factors is given 
by 



H 2 (t)^H 2 (t)-2H 2 (t) 



Sj(t) 

a(t) 



2H (t) 



Si(t) 
a(t)' 



(52) 



where i = 1,2, 3. 

For the volume V(t) of the Universe we obtain 



V(t) = a 3 (t) 



1 



a(t) 



(53) 



and 1/V{t) « l/a 3 {t). 

Hence, in the first order approximation, the gravita- 
tional field equations describing the slightly perturbed 
isotropic flat FLRW Universe become 



3H (t) + 3H 2 (t) - \H (t) + H 2 (t) 

-- Pde + Pdm + 3 (pde +Pdm) 



a(t) a(t) 



+2/3 (p DE + Pde + Pdm + Pdm) ,(54) 



1 d 



a 3 (t) dt 



[a 3 (t)H (t)} 



1 d 
a 3 (t)~dl 



[a 2 (t)H (t)S + (t)} 



+ - 



l_d_ 
?dt 



a 2 (t)Si(t) 



Pde - Pde + Pdm - Pdm 



+0 (pDE + PDE + PDM + PDM) ,(55) 



and 



^[^t )m) ] + ^[At)mm)] 



i__d 

~a?di 



a 2 (t)S 3 (t) 



1 



(pde - Pde + Pdm - Pdm) 



(56) 



where we have denoted 6+(i) = S±(t) + S^i). Equation 
(|46p , giving the time evolution of the volume of the Uni- 
verse, takes the form 

d 2 3 

-li2 a3 ( t ) + 9 l a (*) (P DE ~ P° E +Pdm ~ Pdm)} + 

d 2 3 

-7-^a 2 (t)8(t) + - [a 2 (t) (pde - Pde + Pdm - Pdm)} x 



dt 2 



1 



<*(*) - 5 [ a *®P (Pde + PDE + PDM +Pdm)} = 0. (57) 

The gravitational field equations of the background, 
isotropic geometry, are the standard Friedmann equa- 
tions, determining the scale factor a as 



, a 2 (t) 

\ 2 (t) 



Pde + Pdm, 



a(t) a 2 (t) 

2 ~i7\ "• 77V = ~Pde - Pdm ■ 

a(t) a(t) 



(58) 



(59) 



respectively. With the use of Eqs. (1551) and ([59"f we im- 
mediately obtain 

1 d [a 3 (t)H (t)} I 

a^fj dt " 2 {PDE ~ PDE + PDM ~ Pdm) ' 

(60) 



3H (t) + 3H (t) = -- [pde + Pdm + 3 (pde +Pdm)} 



and 
d 2 a 3 (t) 3 



(61) 



dt 2 



- [a 3 (t) (pde - Pde + Pdm - Pdm)} = 0, 



(62) 
respectively. 

Therefore, the gravitational field equations describ- 
ing the time evolution of a cosmological model with a 
small anisotropy in the z-direction, induced by the non- 
comoving nature of the dark energy and dark matter, 
take the form 



S(t) a(t) 5(t) 
a(t) ~ 'ajtj'^it) 



= ~/3 (pde + Pde + Pdm + Pdm) , 



(63) 



^|[° a (*W*)M*)] + ?|[° 9 ^W 



1 



P (Pde + Pde + Pdm + Pdm) , (64) 



2 d 

a^di 



[a 2 (t)H (t)8i(t)} + 



1 d r 
~^~dl 



a 2 (t)5 3 (t) =0. (65) 



Equations (p5|) - (|55|) represent a system of three dif- 
ferential equations for the three unknown functions 

(S 1 (t),S 3 (t)J(t))- 
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The deceleration parameter q of the non-comoving two 
fluid anisotropic Universe is defined as 



= d 1 

q dt H(t) 



(66) 



and for the case of the small deviation from isotropy it is 
given by 



q = qo + 



1 d 1 



3 dt aH a 



S(t) 






(67) 



where q = d[l/H (t)] /dt- I = - a(t)d(t) / a 2 (t) is the 

deceleration parameter of the isotropic FLRW Universe. 

Equation (|65p can be immediately integrated to give 



2H (t)5 1 (t) + 5 3 (t) = 



C 



(68) 



where C is an arbitrary constant of integration. In order 
to determine the value of C we assume that at t = to the 
Universe was in an approximate isotropic state, so that 
Hi (to) ~ Ho (to), i = 1,2,3, with a small anisotropic 
perturbation Si (to) — Sot, i — 1, 2, 3. Then from Eq. (|4"9")) 
we have Si (to) — Hq (to)Soi, i = 1,2,3. For the initial 
values of S we obtain 8 (t ) = 2<5i (t ) +5% (t ) = S , while 
5 (to) = H (to) So- Hence the value of the integration 
constant C is given by 



C = 



a 2 (to)[2Ho(to)Soi+Ho(to)S 03 ] 



= a 2 (t ) Ho (t ) S f 



(69) 



By multiplying Eq. (164[) by two, and adding the result- 
ing equation to Eq. (1651) . we obtain 



2 d 
a 3 (t) dt 



[a 2 (t)H (t)S(t)] 



a 2 (t)5(t) 



HP- 



DE 



1 d 

+ Z ?~dt 

Pde + Pdm + Pdm)- 



(70) 



By eliminating the term containing the energy densi- 
ties between Eqs. (|6"3"|) and (ITO)) . we obtain the second 
order differential equation governing the behavior of the 
total perturbation S due to the non-comoving motion of 
the dark energy and dark matter as 

2S(t) + 4Ho(t)S(t) + \Ho(t) + 3H 2 (t)] S(t) = 0. (71) 

Equation (|TTj) must be integrated together with the ini- 
tial conditions 



S (t ) = So, 8 (to) = H (to) S . 



(72) 



The behavior of the perturbations to isotropy induced 
by the non-comoving flow of the dark energy and dark 
matter is determined by the Hubble parameter of the 



isotropic FLRW background. Once the function 8(t) is 
known, the time evolution of the perturbations along the 
z— axis is determined by the first order linear differential 
equation 



S 3 (t) - H S 3 (t) 



C_ 

am 



- H S(t), 



with the general solution given by 

C H Q (t) 



S 3 (t) = a(t) { d + 



a 3 (t) a(t) 



S(t) 



dt 



(73) 



(74) 



where C\ is an arbitrary integration constant to be de- 
termined from the initial condition ^3 (to) = Sq 3 - The 
perturbation from isotropy along the x and y axes can 
be obtained from 



1 



S 1 (t) = -{S(t)-S 3 (t)}, 



(75) 



while the anisotropy parameter /3 can be obtained as 

2H Q (t)S(t) + \?,Ho(t)/2 + 5H 2 (t)/2] S(t) 
#*) = -TZ hz ~„ ~Z T ■ (76) 



a (pde + Pde + Pdm + Pdm) 

Hence we have obtained the general solution of the 
perturbed gravitational field equations due to the non- 
comoving motion of the dark energy and dark matter. 

As simple applications of the present formalism we 
will consider two cases, namely, the case of the constant 
Ho, corresponding to an exponentially expanding Uni- 
verse, a(t) = exp (Hot), and of the power law expanding 
isotropic FLRW Universe, with a(t) — t n . Note, however, 
that in order to obtain the expression of the anisotropy 
parameter /3, the energy densities and pressures of the 
dark energy and dark matter must also be known. 



A. Specific solution I: a(t) = exp (Hot) 

If Ho = constant, then the general solution of Eq. ([7TJ) 
is given by 



«*(*) 



-,-Hat 



C 2 cos (V8H t) + C 3 sin ( VsH t 



where C2 and C 3 are arbitrary constants of integration to 
be determined from the initial conditions S (t Q ) = Sq, and 
S (to) = H (t ) So, respectively. For an exponentially ex- 
panding background Universe the magnitude of the total 
anisotropy S decreases exponentially in time. 

By using these expressions of S(t) the perturbations 
Si(t) and S 3 (t) along the x, y and z axes can be explicitly 
computed for these particular cosmological models. For 
the de Sitter-type background expansion, we obtain 



S 3 (t) = C x e Hat + 



- 2g »' {H e H ^ [(C 3 - V2C 2 ) sin (2V2H t) + (C 2 + y/2C 3 ) cos (2y/2H t)] - 2C) 



and 



Si(t) 



-2H t |2 ( c _ 3CiHof^*) + g e got [(5C 2 - y/2C 3 ) sin (2V2H t) + (5C 2 + V2C 3 ) cos (2>/2ff t)] } 

12#o 

I 



(78) 



(79) 



respectively. The time evolution of the anisotropic per- 
turbation along the 2— axis is proportional to the scale 
factor of the isotropic background, and thus for a de 
Sitter-type Universe the anisotropy in the z direction 
will increase in time, making the Universe more and 
more anisotropic. On the other hand, since in the large 
time limit 8 tends to zero, lim^oo 8 — 0, from Eq. ([75)1 
it follows that lim^oo 8\(t) = — (1/2) lim t _j. oc 8 3 (t) = 
— (1/2) limt^oo C\ exp(H t). However, in these large 
time limits the formalism developed in the present paper 
may not be applicable, since all the presented results have 
been derived under the assumption Si <C a, i = 1,2,3. 



B. Specific solution II: power law expansion 

For a power law expansion, with the scale factor of the 
form a(t) = t n , Eq. ([TTj) becomes 



•j in ■ n (3n — 1) .,,. 
25+— 8+ — ^ >-5(t) = 0, 

with the general solution given by 

S(t) = C±t ri + C 5 t r2 , 
where 



n,2 = 



1 - 2n ± V2V2 



(80) 



(81) 



(82) 



In the case of the power law background expansion we 
obtain 



V. OBSERVATIONAL IMPRINTS ON THE CMB 

AND SUPERNOVAE LUMINOSITIES OF THE 

NON-COMOVING DARK MATTER-DARK 

ENERGY MODEL 

In the present Section, we briefly consider the ob- 
servational implications of the non-comoving dark en- 
ergy and dark matter models. The anisotropic expan- 
sion of the Universe would leave its imprints on both 
CMB and supernovae date. The amount of possible 
anisotropy existing in the Universe can be constrained 
by using the CMB and SNIa observations. Usually the 
present CMB data constraints can be much tighter than 
the constraints from the anisotropies in the luminosity- 
distance-redshift relationship of SNIa [23|]. However, 
there are anisotropic models which can avoid this bound 
completely .23]. Within these bounds, the existence of 
anisotropy in the Universe can provide a potential expla- 
nation of various anomalies in the cosmological observa- 
tions, especially in the CMB at the largest angles. 



A. Effects on CMB on the non-comoving motion of 
dark energy and dark matter 

A possible anisotropic expansion of the Universe will 
leave its imprint in the CMB. To calculate the CMB spec- 
trum in a background geometry with a Bianchi type I 
metric given by Eq. (|29[) we follow the approach intro- 
duced in [23]. Photons emitted by astrophysical sources 
at cosmological distances move in the Bianchi type I Uni- 
verse along the geodesic lines, given by 



5 3 (t)=dt n + t 



C 4 H t r i C 5 H t r2 



ct~ 2n 



1 — 3n n — r\ — 1 n — ?*; 



1 



while for <5i(t), we have 

2*i (t) = df' 1 ~ C x t n + +C 5 t r * + C 
H a t ri+1 



t l 



+C 4 - 



-n + n + 1 



C s 



3n-l 

H t r2+1 



-n + r 2 + 1 



(83) 



(84) 



The arbitrary integration constants C4 and C5 can be 
determined again from the initial conditions 6 (to) = 5q, 
and 6 (to) = Ho (to) So, respectively. 



~d\ 



iyv=o, 



(85) 



where A is an affine parameter, and the Christoffel sym- 
bols r^ s associated to the metric Eq. ([2Uf are given by 

Mi = a i^i' an d To* = Hi, i = 1,2,3 (no summation 
upon i in the Christoffel symbols). 

The four- velocity u M = dx^ /dX of the photons satisfy 



the normalization condition 



= 0, giving (uqY 



af (u l ) . Consider two photons, emitted at two time in- 
terval (to = t e , t\ = t e + r), where t « t c . By taking the 
difference of the two normalization conditions for pho- 
tons in the first order of approximation in r we obtain 
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, dr 



^^,,( M fr(A) + 0(r 2 ). (86) 



d\ 



;=i 



By introducing the redshift z, defined as 1 + z (A e ) = 
r (A r ) /r (A e ), where r (A r ) is the time difference of the 
received signals, with the use of Eq. (|56| we immediately 
obtain 



d_ 
d\ 



1 3 



(87) 



The photon velocity components u l , i = 1, 2, 3, can be 
obtained from the photon geodesic equation of motion, 



aA ai 



giving 



«'(*) 






* = 1,2,3, 



* = 1,2,3, 



(89) 



where uoi, i = 1, 2,3 are arbitrary constants of integra- 
tion. In order to determine the integration constants, we 
adopt the convention according to which the present day 
t = to value of the scale factors is ai (to) = 1. We de- 
note the present day values of the photon four-velocity 
as as u l (to) — u 1 , i — 1,2,3. Since one can always repa- 
rameterize the affine parameter A without modifying the 
physical results, we can normalize the present day pho- 
ton four- velocities so that Yli=i (v 1 ) — 1- Hence one 
may interpret the unit vector u in terms of the angles 
(u x ,u y ,u z ) = (sin 9 cos 0, sin 9 sin <f>, cos #), representing 
the present day arrival angles of the photons to the ob- 
server. Therefore, by substituting the photon velocities 
into the redshift equation, we obtain [23| 



respectively. 

The temperature distribution is determined by the re- 
lation 



T 



(*)- T 



T, 



+ z [u 



(94) 



where the last scattering temperature T* does not depend 
on the direction. However, due the anisotropic nature of 
the geometry of the Universe, photons coming from dif- 
ferent directions will be redshifted by different amounts. 
The spatial average T of the temperature field is obtained 

as AttT = J T (u) d£l~, while the anisotropies in the tem- 
perature field are given by 



ST 



(») 



T 



(*) 



T 



= 1 - 



T* 



T v /l + 2[J 1 (i )-<J 3 (to)]cos2e 



(95) 



The multipole spectrum C; can be described by the 
coefficients in the spherical expansion of the temperature 
anisotropy field. The quadrupole Qi can be obtained as 

231 ' ' 



C-2 



5\/3 

4 1 



-z~y 



5%/3 



dVZa 



Mt) - 5 3 (t)} 



(96) 



The nominal value of the best fit quadrupole measured 
by Planck is [2(2+l)/27r]C 2 = 299.495x 1(T 12 @. There- 
fore the recent observational results can be used to con- 
strain the parameters of the non-comoving dark energy- 
dark matter model. 



1 + z 



(*) 



3 -2 
Uf 



^s« 



which can be rewritten as 

1 
«i 



' '< z[u) = —\ll + ulel + u%el, 



where the eccentricities e 2 and e 2 z are defined as 



a 2 



^'=[21 | I. 

a-3 



(90) 



(91) 



(92) 



Since in the non-comoving dark energy-dark matter 
model we have considered we have assumed that a\ =02, 



we obtain e = 0, and 



Mt)-s 3 (t)}, 



(93) 



B. Constraints on non-comoving dark matter and 
dark energy from supernova observations 

Anisotropy in the expansion rates of the Universe 
would, in principle, determine rotationally non-invariant 
luminosity distance-redshift relationships for the la type 
supernovae. The luminosity-redshift relationship of the 
SNIa can be used to probe the possible anisotropies in the 
expansion history of the Universe [2J, [2J] . The spatial 
geodesic Eq. ([55)1 shows that the direction of a photon u 
that is coming towards us from an astrophysical source is 
constant. This follows from the orthogonality property 
of the Bianchi I models. The null geodesic equation can 
be rewritten as 



j,2 1 -2. 2 1 -2 2 1 ~2 2\ 1 
at = \u x a 1 + u y a 2 + u z a 3 ) dr 



(97) 



y 2 + z 2 . The conformal distance can be 



where r 2 = x 

obtained as J dr. The luminosity distance at the redshift 
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z in the direction u is thus given by [23!, \M 
'*<*> alt 



d.L (z,u) = (1 + z) 



o ^u 2 x a\ + u 2 y a\ + u\a\ 



(98) 



Therefore high precision SNIa data can be used to con- 
strain the non-comoving dark energy-dark matter mod- 
els. 



VI. DISCUSSIONS AND FINAL REMARKS 

In the present paper, we have considered the theoreti- 
cal possibility that our Universe may be modeled as a 
mixture of two non-interacting perfect fluids, namely, 
dark energy and dark matter, respectively, with differ- 
ent four-velocities. This model is formally equivalent to 
a Universe filled with a single anisotropic fluid, with dif- 
ferent pressure components along the coordinate axes. 
The gravitational equations describing the non-comoving 
dark energy and dark matter mixture were explicitly ob- 
tained for a Bianchi type I geometry. By assuming that 
the difference in the four-velocities of the dark energy 
and dark matter is small, one can consider that the de- 
viations from isotropy are also small, representing just 
a small perturbation of the homogeneous and isotropic 
background metric. In particular, the perturbation equa- 
tions have been explicitly obtained and solved, for the 
case of the Bianchi type I geometry. Two simple cases, 
corresponding to an exponential and power law expan- 
sion of the background geometry have also been analyzed 
in detail. 

Therefore, if there is a slight difference between the 
four-velocities of the dark energy and dark matter, the 
Universe would acquire some anisotropic characteristics, 
and its geometry will deviate from the standard FLRW 
one. In fact, the recent Planck results show that the pres- 
ence of an intrinsic large scale anisotropy in the Universe 
cannot be excluded a priori, so that the model presented 
in this work can be considered as a plausible and vi- 
able working hypothesis. However, in the present paper 
we have considered that the deviations from isotropy are 
small, and we have restricted our analysis to the linear 
perturbation regime. If the amplitude of the perturba- 
tions does not satisfy the condition Si <§C a, i — 1,2,3, 
then non-linear effects should also be taken into account, 
and incorporated into the description of the anisotropic 
effects generated by the non-comoving dark energy and 
dark matter. 

We have also pointed out some possibilities of observa- 
tionally testing the model, by using the recently released 
Planck data [l(, as well as the type la supernovae data 
[3| • In the case of the CMB we have explicitly obtained 
the distribution of the temperature anisotropy and of the 
quadrupole C*2 as functions of the deviations Si and S3 
from the isotropic flat geometry. Once the form of the 
background FLRW geometry, more exactly the form of 
the scale factor a(t) is known, one can use the obtained 



relations to put tight constraints on the parameters of 
the non-comoving dark energy-dark matter model. Sim- 
ilar constraints can be obtained from the study of the 
anisotropic luminosity distance by using the type la su- 
pernovae observational data. 

In this context, we emphasize that a major challenge 
for cosmology is to elucidate the nature of dark energy, 
which drives the accelerated expansion of the Universe. 
The simplest explanation for dark energy is the cosmo- 
logical constant, where the latter has an equation of state 
w = pde/pde = — 1- Scalar fields, alternatives to the 
cosmological constant, have an equation of state with 
w > —1. However, the CMB data alone does not strongly 
constrain w, due to the two-dimensional geometric de- 
generacy in the theoretical models. The recent results 
of Planck, combined with the supernova data, provide 
w = -1.09 ± 0.17, or w = — 1.13±g;i|, generally in a 
good agreement with the cosmological constant hypoth- 
esis. However, at the present level of observational pre- 
cision one cannot rule out the existence of a dark energy 
component of the Universe, distinct from a cosmological 
constant. For the case of a "pure" cosmological constant, 
the formalism developed in the present paper cannot be 
applied, since there is no flow associated to A. 

In the case of the scalar field <f>, with energy-momentum 
tensor 






V M </>V„</> - 5 „„V A </>V A </>/2 - U(0) 5ai1 ,, (99) 



where V(4>) is the scalar field potential, a perfect fluid de- 
scription corresponding to an energy-momentum tensor 
Tde ilv = (P0 +P4>)u (t> M U0 v - p^g^y can be obtained, if 
V A </> is timelike, by introducing the effective four-velocity 
of the scalar field fluid as 12511 



Ufa 



y/\V x <F7x<f>\ ' 



(100) 



where V A </>V A <^> ^ 0, and UfaU§ = sign (V A 0V A </>). The 
effective fluid pressure and energy density associated to 
the scalar field energy momentum tensor can be obtained 
as P<t> = TJ$v,$v% and P4> = T^h^/3, where hf" = 
g^ v — u'iu'^, and are given by |25| 

P4> = [V A 0V A </>/2 - V{<f>)} sign (V A 0V A 0) , (101) 



and 



= [-l + sign(V A (/»VA0)/2]V A ( /»V; 
-[4 + sign(V A ^VA^)]Vfa), 



(102) 



respectively. Therefore if the four-velocities of the dark 
matter fields, u^, M and of the scalar field vP satisfy the 
condition 



L DM 



*<, 



(103) 



the mixture of the non-comoving scalar field (dark en- 
ergy) and dark matter can be described by a single 
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anisotropic fluid. Such a situation may have occurred 
during the inflationary, or post-inflationary (reheating) 
era of the evolution of the Universe. 

The present model can be easily generalized for the 
case of a three component mixture, or for cases involving 
the presence of condensates and multi-component scalar 
fields. These possibilities, their impact on the cosmologi- 
cal evolution, as well as some observational implications, 
will be considered in future studies. 
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